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Abstract: We study non-gaussianity effects, using the 5N formalism, in a multi-field in- 
flationary model consisting of Kahler moduli derived from type IIB string compactification 
in the large volume limit. The analytical work in this paper mostly follows the separable 
potential method developed by Vernizzi and Wands. The numerical analysis is then used 
in computing non-gaussianity beyond slow-roll regime. The possibility of the curvaton 
scenario is also discussed. We give the condition for the existence of the curvaton and 
calculate the non-guassianity generated by the curvaton decay in the large volume limit. 
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1. Introduction 

Since it was first proposed in the early 1980s, inflation has provided an important insight 
into the understanding of the very early universe [1, 2, 3]. It also affords the possibility 
of probing the fundamental theories that provide a microscopic explanation of inflation 
such as string theory. (For a review, see [4, 5].) In recent years, there have been some 
promising developments of models derived or inspired by string theory, such as type II 
flux compactifications [6, 7], where scalar fields associated to the shape and size of the 
internal space, or to the positions of branes, serve as candidate inflaton fields [8]. There 
now exists many inflationary models derived from string compactifications which correctly 
can account for the observational results from WMAP, such as the spectral index, n^, and 
the ratio of tensor to scalar perturbations, see e.g. for a review [9]. 

Due to the nature of cosmological models arising from string compactifications, these 
inflationary models have several interesting features: (i) in general, the scalar potential 
is a highly non-trivial function that depends on many scalar fields and (ii) the fields are 
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not canonically normalized, i.e., the metric is in general neither diagonal nor field inde- 
pendent. These features of string theory lead to a highly coupled dynamical situation, in 
which, in principle, the motion of any one field impacts the evolution of the other fields. 
Therefore, we need to understand how each modulus and its perturbations evolve during 
and after inflation, including both light and heavy fields. In this paper, we extend our ear- 
lier study [10] to include the non-gaussian fiuctuations in a multi-field infiationary model 
arising from string theory. 

Gaussian fluctuations are described by the two-point function and the corresponding 
power spectrum. The observed primordial perturbation to date is gaussian to a good ac- 
curacy in agreement with the predictions of inflation. However, it is possible that future 
experiments may observe non-gaussianity in primordial perturbation. It is therefore im- 
portant to study the non-linear effects in inflationary models that can can give rise to such 
non-gaussian fluctuations. 

In single field infiationary case, the non-linear parameters /tvl and qnl, which char- 
acterize the size of non-gaussianity, can be calculated in terms of the slow-roll parame- 
ters [11] [12]. The result is generally small, of the order of slow-roll parameters. In multi- 
field models, there are usually both heavy fields and light fields. The light fields are believed 
to drive inflation and heavy fields are frozen during inflation. It is often assumed that the 
dominate contribution to non-gaussianity comes from the inflaton. However, for an arbi- 
trary scalar potential, it is not clear that whether the other (non-inflaton) fields have any 
sizable contribution to non-gaussinity. We would like to address this issue, at least for the 
model studied in this paper. 

In the infiationary scenario in which the primordial curvature perturbation originates 
from the inflaton, other light flelds only play a role in assisting with stabilizing the poten- 
tial. In the curvaton scenario, however, the curvaton (light, non-infiaton field) may have 
significant contribution to the primordial perturbation if its energy density grows large 
enough at a later time after inflation but before it decays into radiation. We explore the 
possibility of a curvaton scenario and compute the amount of non-gaussianity generated 
by the curvaton. 

We focus on a string inspired model based on the large volume scenario [13] [14] [15] . We 
study different conflgurations of the model with different numbers of scalar flelds (moduli) 
and with various values of the volume of the Calabi-Yau to provide hints of the microscopic 
physics by connecting non-gaussianity (if observable) with the model parameters (moduli, 
the volume, etc). 

The outline of the paper is as follows. We first review the non-gaussian perturbations 
and the 6N formalism in Section 2. In Section 3, we introduce the scalar potential arising in 
the large volume scenario of type IIB string compactifications. We then apply the separable 
potential method to the above multi- field infiationary model in Section 4. In Section 5, a 
numerical analysis is carried out that extends the previous analytical study beyond slow- 
roll. Comparing the two methods, we find a good agreement in the regions where they 
overlap. Finally, in Section 6 we study under what conditions a curvaton may exist after 
the end of infiation in this type of model derived from string theory, and calculate the 
contributions to /nl- 
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2. Non-gaussian perturbations and the 5N formalism 

In the following, we briefly review some general facts about non-gaussian perturbations 
and the 5N formalism which is a powerful tool to calculate non-gaussian effects. (For more 
detailed discussions, see eg [11] [16].) 

Let us define the e-folds time ^ 

N = j^'^ Hdt (2.1) 

where t* is usually chosen to be some time during inflation (the initial flat slice) and tc is 
some epoch later with constant curvature perturbation (the final slice of uniform density). 
In multi-field theory, the primordial curvature perturbaiton reads [16] [17] [18] 

C(x) = 6N = NM + ^NijS^'6ct>^ + ^N,jkS(p'6(l^6^'' + ... (2.2) 

where Ni,Nij, ... are derivatives of the e-folds with respect to the fields 6(p'^ are evaluated 
on the initial (flat) slice while the derivatives of N are evaluated on the unperturbed 
trajectory with respect to the unperturbed fields at Hubble crossing[19]. In the later 
sections, we will use to denote Hubble crossing. 

The three-point correlation function is given in terms of the bispectrum B(^{ki, k2, k^) 

(Ck.Ck.Cka) = (2vr)353(ki + k2 + ks)B^{ki,k2, ks) (2.3) 

where Ck; are the Fourier coefficients of C(x). If slow-roll is satisfied, the bispectrum is 
completely specified by the non-linear parameter fNLiki, k2, k^) [20] [22]: 

6 

B(;{ki,k2, ks) = -fNL{h,k2, ks)[P(;{ki)Pi^{k2) + cyclic permutations] (2.4) 

where P({ki) is the power spectrum. Using the 5N formulation, the non-linear parameter 
fNL is given by [23] [24] [25] 



-^fNL{h,k2,h) = ,3 ^ ,3 ^ ,3 4,4 P2 = ^^^2 + 0{r) (2.5) 



where r is the tensor to scalar ratio. The correction 0{r) is a fcj-dependent geometric 
term [25] [26]. In standard slow-roll inflation, 

r~16e, e~--^<l (2.6) 

Thus 0{r) is much less than unity due to slow-roll condition and observation constraints 
on r [27] [28]. From now on we focus on the first term in (2.5) and redefine the momentum- 
independent non-linear parameter 



fNL = ^ ^';^ (2.7) 



The four-point function has the form 

(CkiCk.Ck3Ck4) = (27r)353(ki +k2 +k3)T^(ki,k2,k3,k4) (2. 
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Neglecting corrections of the order of the slow-roh parameters, the trispectrum reads [11] [20] [21] 



T^(ki,k2,k3,k4) =TNL[P({ki)P(ik2)P(:ikiA) + H permutations] 



54 

+ :^9NL[P({k2)P(;{k3)P(:{ki) + 3 permutations] 



25 



where ki4 = |ki — k4|. Note that, unhke the bispectrum, the trispectrum depends on the 
directions of /cj's. The parameters, when corrections of slow-roll order are neglected, are 
given by [18] 

'^NL = (2.9) 



25 NijkN'N^N^ 
^^^ = 54 {N.N^Y (^-^^^ 

Applying the Cauchy-Schwarz inequality to (2.7) and (2.9), we get the following relation [29] 

TNL > f^flL (2.11) 

3. The Scalar Potential and Kahler Moduli Stabilization 

In what follows we focus on a particularly inflationary model derived from string theory 
consisting of multiple Kahler moduli, in the large volume limit (also known as the Large 
Volume Scenario) [13] [14]. 

Supergravity in a four dimensional theory with N = 1 supersymmetry is completely 
specified by the real Kahler potential JC{ip^, (^*) and the holomorphic superpotential W{ip^). 
Focusing on the dynamics of the scalar fields relevant for inflation, the supergravity action 
is (we will work in the Einstein frame, and in units where Mp = 1) 



Sx=i = d x^f^ 



(3.1) 



The scalar potential depends on the superpotential W, the Kahler potential K as well as 
the Kahler metric G^j, 

V{cP^, = (g'WiWD-W - 3 w) + Kpua (3.2) 
DiW = diW + WdiK (3.3) 
= did-^K (3.4) 

The derivatives di and di differentiate with respect to the ipi and (pi dependence, respec- 
tively. By expanding the complex fields in terms of their real and imaginary part, we 
can relate the supergravity action above, (3.1), to the action, (3.13). The term VupUft will 
include the effects of supersymmetry breaking arising from other sectors of the theory. 

We will demonstrate our methods in the context of Type IIB string theory compactified 
to four dimensions on a Calabi-Yau orientifold because the scalar potential in this case is 
well-understood and realistic four-dimensional models can be constructed [6, 7, 13, 14, 30, 
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31]. After including the leading perturbative and non-perturbative corrections of string 
theory, the Kahler potential and superpotential are given by 



3 



/C = -2 In I V + j - hi(-f(r - f )) - hi j SI A , 

|^^G3A0 + J^Aie-'^'^') (3.5) 



w J 



3 (3.6) 



Here Qs is the string couphng, Ig is the string length, Q is the holomorphic three-form on 
the Calabi-Yau manifold M, G3 is the background field (flux) that is chosen to thread 
3-cycles in M and 

C(3)x(M ) 
^ 2(27r) 

where x is Euler number of M. The axion-dilaton field is t = Co + i e"*^, and the integrals 
involving Q are implicitly functions of the complex structure moduli. The fields Tj = 
Tj + ibi are the complexified Kahler moduli where Tj is a 4-cycle volume (of the divisor 
Di G Hi{M^'L)) and bi is its axionic partner arising ultimately from the 4-form field. 
Here aj = 271 /Ni for some integer Ni, for each field, that is determined by the dynamical 
origin of the exponentials in the superpotential [Ni = 1 for brane instanton contributions, 
A^i > 1 for gaugino condensates). Finally, V is the dimensionless classical volume of the 
compactification manifold M (in Einstein frame, but measured in units of the string length). 
In terms of the Kahler class J = (by Poincare duality Di € H'^{M,1i)), with the 

measuring the areas of 2-cycles, Cj, 



V = / = -KijktHH'' , (3.7) 
Im 6 



where Kij^ are the intersection numbers of the manifold. V should be understood as an 
implicit function of the complexified 4-cycle moduli via the relation 

n = duV = ^K.jktH'' . (3.8) 

There are additional perturbative corrections to /C in (3.5), but we have kept the terms 
that give the leading contributions to the scalar potential in the large V limit of interest 
to us [32]. In particular, expanding fC to linear order in ^ gives a consistent approximation 
in inverse powers of V. We have also assumed that all of the Kahler moduli Tj appear in 
the superpotential (see [30] for examples) and that we use a basis of 4-cycles such that 
the exponential terms in W take the form exp{—aiTi). As these exponentials arise from 
an instanton expansion, in order to only keep the first term as we have done, the 4-cycle 
volumes must be sufficiently large to ensure that ajTj ^ 1. 

Finally, the form of the term KipUft in (3.2) depends on the kind of supersymmetry 
breaking effects that arise from other sectors of the theory. We take 

Kxpiift = ^ (3.9) 
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which will describe the energy of a space-filling antibrane [7], fluxes of gauge fields living 
on D7-branes [33], or the F-term due to a non-supersymmetric solution for the complex 
structure/axion-dilaton moduli [34]. 

It was shown in [7] that a generic choice of background fields G3 causes all the complex 
structure moduli and the axion-dilaton to acquire string scale masses without breaking su- 
persymmetry. They are then decoupled from the low-energy theory and their contributions 
to fC and W are constants for our purposes^: 

IC = -2ln(v + i] - In ( -] + /Co, 
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gs 



W = -^iWo + Y, A^e-""^^^ J , (3.10) 



where /Co (Wq) is the complex structure Kahler potential (superpotential), evaluated at 
the locations where the complex structure moduli have been fixed. It was shown in [13] 
that, when the Euler number, x < 0, for generic values of Wq (and hence of the background 
fluxes 6*3 ), the scalar potential for the Kahler moduli has a minimum where the volume V 
of the Calabi-Yau manifold M is very large - the associated energy scale is a few orders 
of magnitude lower than the GUT scale. Furthermore, in these Large Volume Scenarios 
there is a natural hierarchy - one of the Kahler moduli is much larger than the others and 
dominates the volume of the manifold, 

Ti > T2,r3,r4 • • • (3.11) 

which we will use to simplify the effective potential. For our purposes these models are 
also attractive because the scalar potential admits an expansion in inverse powers of the 
large volume V. This will allow us to carry out analytical calculations of inflation arising 
from Kahler moduli rolling towards the large volume minimum of the potential. 

For transparency of the equations, we will assume that the intersection numbers kijk 
are such that in the basis of 4-cycles, r^, the volume takes the diagonal form [14] 

V = a(Tii -^AiT^i) = -a^A,r,i (3.12) 

i=2 1=1 

where Ai = — 1, and Aj,z > 2 are usually positive. 

With the volume taking the above form we can explicitly compute the metric on the 
moduli space, Qf^ = didjlC, which is needed both for the kinetic energy ^h^^fifj, where 
hij = "iQfj, and for the scalar potential, F, appearing in the four dimensional action, 

s = j d^x^[^ - h^^d^nd^^T, - y(T,)] (3.13) 



^In the case of the F-term breaking due to the complex structure/axion-dilaton moduli [34], the con- 
tribution of the complex structure and axion-dilaton moduli to the scalar potential does depend on the 
volume (3.9). 
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By expanding in inverse powers of V, keeping terms to 0{V ), we obtain 

/i - . = ^ ; , 6'^ + . (3.14) 

With the axions minimized in the potential, the effective potential then becomes [15] 

i=2 i=2 

where we have assumed that /Co can be chosen such that the overall scale of the potential 
is simplified, i.e., overall factors of gg and 27r are not present. Here we have expanded V 
to 0{V~^) to include the leading a'-corrections, ^^^3 , as well as the uplift term, The 
parameters in the potential can be chosen and tuned under certain constraints [35, 36, 37, 
38]. 

To canonically normalize the metric, we can apply the following field transforma- 
tions [10] 



3Ai(l + 3Ai) 

1 log(Ti) (3.16) 



^' = \h^nK i>2 (3.17) 
V 3ri2 

by keeping terms to leading order in the expansion of inverse powers of the volume, and 
using that in the large volume scenario 

V^an^ . (3.18) 

This field redefinition results in a good approximation to the canonical metric. After 
canonical normalization, the equations of motion read 

dV 

cj)' + 3H(l)' + —^=0, i = l,...,n. (3.19) 

To get successful inflation, we take the following steps (for more detail, see [10]). First, 
find the global minimum of the potential by 

dV 

— = 0, i = l,...,n. (3.20) 

Written explicitly, 

dV 

V=?^i^yn.-, i = 2 „. (3.22) 

where (3.21) is obtained by 



dV dV dV , , 
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and we apply (3.18) 

dV 

7^ = 0, i>2 (3.24) 

OTi 

to get (3.22). To have a small but positive cosmological constant, we also require that 

Knin > (3.25) 

In practice we want ajTj S> 1 so that all higher order non-perturbative corrections of the 
form ^-^"'i'^i ^ with integer m > 1, in the scalar potential are negligible and the effective 
potential (3.15) becomes a valid approximation. We find that the global minimum of the 
potential only exist with the parameters lying in certain regions of the parameter space. 
Let the global minimum be 

^^iii) ''"imin- (3.26) 

We have 

Vmin = (3.27) 

where 

n 



i=2 al 

and 



= (3.29) 

Then, for inflation to start, we displace the flelds away from the global minimum along 
the flat direction of the potential, and find the corresponding local minimum. Denote the 
values of the fields at the local minimum by 

VinUTiim- (3.30) 

which are the initial conditions of the model. By solving the equation of motion (3.19), we 
should get successful infiation in which the fields evolve slowly toward the global minimum 

(^imin) • 



4. The Inflaton Scenario 

4.1 Separable Potential Method 

In order to use the 5A^-formalism, we need to calculate the derivatives of the number of 
e-folds, N , with respect to the fields. For an arbitrarily shaped potential, this can be done 
using numerical method which will be discussed in the next section, while the analytic 
treatment of non-gaussianity is known for being difficult. If the potential satisfies certain 
criteria ([25], [39], for example), the e-folds can be obtained by analytical integration . In 
the model introduced in Section 3, it has been shown [40] [10] that the non-infiaton (or 
heavy) fields are frozen before the end of infiation and only the light fields evolve during 
inflation. Thus, before the end of inflation, we can ignore the dependence of the potential 
on the heavy flelds and treat the potential as a function of the light flelds only. 
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Let us assume that we have a volume modulus, ti, a heavy modulus, T2, and two light 
moduli, T3 and T4 (as in [10]), corresponding to the canonically normalized fields (p^ and 
04, (3.16)-(3.17). Furthermore, if we suppose that ri and T2 are frozen during inflation, 
the potential (3.15) can be seperated as two functions each depending only on one of the 
light fields {(j)3,(j)4): 

V = UicPs) + W{(l)4) (4.1) 

The terms which contain the frozen fields have been absorbed into U{(j)-i) and W{(j)4). 

Next we will follow the separable potential method developed by Vernizzi and Wands [25] 
to calculate the derivatives of the e-folds. In the canonical frame, assuming the background 
fields only have time dependence, the background equations of motion read 

+ 2,H(i)^ + y." = 0, a = 3, 4. (4.2) 

where 



W\ a = 4. 



And the Einstein field equations are^ 



^H^ = -g^'d^n^a + V (4.4) 



2 

Using the slow-roll approximation, we have 



H = -Ig^'d^rOucl^a (4.5) 



o V 

H^^- (4.6) 

1 3(p3 3^4 

H ^ U' ^ W 



By integrating (4.7) we then get 



The number of e-folds becomes 



tt2 

N= I —dt 



-dt+ / dt 



(4.7) 



jw-jw-^ (4.8) 



3H 3H 

Here * and c denote the time at Hubble exit and some time after Hubble exit, respectively. 
We usually choose the latter as the time, tc, for some constant Hubble parameter (uniform 



Again we set the planck mass Mi — 1. 
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energy density hypersuface) , He- Then we let vary and compute the derivative of 

with respect to the initial fields at t*. The results, derived by Vernizizi and David [25], are, 
Ni, = dN/dct>%=u, 

N33* = 1 - f^^^^ + 4.12 

_ m* W^-Z^ 1 dZc 

iV44^, - 1 - — -—==——, (4.1Jj 

1 1 9Zc ,,,,, 



where 



and 



l/KaV 7 ^£3 - ^€4 



v^^r— = -vei^^r— = -v2K— 2 ^ 1 ] 4.16 

Note that this method rehes on the slow-roll approximation (4.6). It requires that 
the final slice at tc must lie within the slow-roll regime. To calculate the amount of non- 
gaussianity generated near and after the end of inflation, one needs to find an alternative 
method valid beyond slow-roll. For example, the authors of [39] proposed an analytic 
method, valid for certain classes of inflation models with separable Hubble functions, which 
can be used to study non-gaussianity after inflation ends. Although their analysis applies to 
certain types of potentials with exponential terms, the detailed conditions are not satisfied 
for the scalar potential (3.15). In section 5, we present a numerical analysis valid beyond 
the slow- roll regime. 

4.2 Estimate of /tvl 

Two light fields case: In the model discussed in Section 3, the two light fields (i;^3, 
serve as candidates for inflaton. We know make a rough estimate of f^i^ in this case. Let us 
assume that is the assisting field and thus lighter than the inflaton <^ . During slow-roll 
we would then expect that 

€4 < £3 < l,lr/4l < ]%! < 1,|??»| >ei, i = 3,4 (4.17) 

since it follows that \y^ii\ ^ Vi from the flatness of the potential (3.15) and the fact that 
both fields are light, \r\i\ = J-^J ^ 1. Let 

\W{4>%=v\V{<l>'-)\. P>0, (4.18) 
£4 = ses, < s < 1, (4.19) 
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and 



q\m\, o<q<i. 



(4.20) 



Although, in general, p, s, q are functions of time, they are not expected to change dramat- 
ically in the slow-roll regime. In particular, we can treat p as a constant due to the flatness 
of the potential. Furthermore, while the actual value of p depends on the specific values 
chosen for the parameters in W{(j)'^) and U{(p^), respectively, one finds that p ~ 0(1) in a 
typical model. Using the separable potential method introduced in the previous section, 
we get 

P 



3* 



1 

\/2e^l + Sc' 

P m* 



4* 



3* 



1 



1 



1 + Sc 2e3* 
ScQ* P m* 



+ 



Sc{Sc + Qc) mc 
£3* 



{l + Sc 



+ 



Sc Sc{sc + qc) mc 



iV34* 



1 + Sc 2e3^ (1 + Sc 

, 1 ScjSc + Qc) mc 

y/s^ (1 + ScY £3* 



£3* 



(4.21) 
(4.22) 
(4.23) 
(4.24) 



The non-linear parameter /atl, defined by (2.7) and evaluated at for a fixed tc, then 
becomes 

Inl ~ XT]?,* + ye3* (4.25) 



where the coefficients are 



-(1 + — 

2 



Sc{,Sc + qc 

p2(l + Sc 



■(1 



Sc-^2mc 
'5* ^3* 



l + Sc 



{l + q* 



2 {l + Scf 

P 1 + ^ 



both of which are of O(^) because of the slowly changing s and q. 

Mixed case: Consider a model in which there is one heavy fields, (^2, and one light field 
03 (inflaton) in addition to the volume field, As usual, 0i is frozen during inflation. 
But we drop the constraint that the heavy field ^2 has been stabilized. We want to see 
how much the heavy field contributes to the non-gaussianity. The potential (3.15) can be 
separated as 

V = i/(<^2) + ^^(03) 



where 



First, since (j)2 is heavy and 

r/2 = 



%2^i2)V^^_2a,. 

3VA2a 



8(03^3)2 



3VA3a 
!»3 is light, 



"^3 g-2a3T3 



402^2^/0 7"2 

4a3A3Wor3 



-a3T"3 



V 



o{i),m 



K33 

V 



<C1; €2 <C €3 



(4.26) 

(4.27) 
(4.28) 

(4.29) 
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Because the inflaton (j)^ is usually displaced far away along the flat direction of the potential, 
we would expect that 

WicPs) < U{(b2) (4.30) 
In addition, we assume that during slow-roll 

ec~e*, r]c^'n* (4-31) 
We then have from (4.10) and (4.11) 

iVa. ~ iVg* ~ (4.32) 



^3* V ^3* 

Since (p2 is heavy and 03 is light it follows from (4.29) that 

A^3* > (4.33) 

and (4.16) can be written 



dz, yejcm^^ dz^^y^2cm^ ^^^^^ 



092* esc V £2* Ctps* esc v ^3* 
Using these results we can estimate the expressions for Nij, and (4.12-4.14) becomes 

iV22*~0(l), iV33~- — , iV23*~— (4.35) 

£3* e3cv^2*e3* 



from which it follows that 



A^33*~ — J— A^23* (4.36) 

mc V «^2c 



By assumption, (4.29), the slow-roll factors satisfy 

^»1, ^»1, — ~0(1), (4.37) 

e2c V ^2* e2c ??2c 

Using (4.37) in (4.36) we then arrive at the following hierarchy among the Nij^ 

N33* > iV23* > iV22* (4.38) 

Therefore, 

„ 5 J:lJ=2N^N,N,^ 5iV33 . ^ ^ _ 

= 6 iEiN^- " 6 « ^ 

with the dominant contribution coming from the light field ^3. 

In summary, both examples discussed in this section yield /atl ~ 0(r/) <^ 1 (where r] 
is the slow- roll parameters for the inflaton). Since \r}\ >> e for the inflaton, the dominant 
contribution to /tvl is indeed the momentum independent term in (2.5). The result is 
similar to that of the standard slow-roll inflation, see for example [11]. 



- 12 - 



4.3 Explicit Setups 



We now perform an explicit calculation of /nl given by (2.7). To evaluate Ni and Nij 
in (4.10-4.16), we need to compute Z at fixed tc- This is done by solving the equations of 
motion for the background fields numerically. Below we present to examples corresponding 
to two different values for the volume of the Calabi-Yau. 

Example 1 Let us choose the parameters in the potential (3.15) as 



1 



2tt 



27r 



2tt 



a 



:,«4 



A2 = 0.2,^3 = 0.001,^4 = 0.001. 



100 100 

Ai = -1, A2 = 0.1, A3 = 0.005, A4 = 0.005, W = 500, ^ = 40, 7 = 9.75 x 10"^ 
Choose the initial conditions to be 

ri(0) = 76212.1, r2(0) = 246.99, T3(0) = 472.42, r3(0) = 491.54, 

fi(0) = f2(0) = 0,7^3(0) = -1.72 X 10^^^f4(0) = -1.5 X 10"^^ 

The volume in this setup is V ~ 10^ which is within a reasonable range 10^-10^ [41]. This 
will give 60 e-folds before the end of inflation. The nonlinear coefficients are ^ 





N=20 


A^=30 


iV=40 


A^=50 


iV=55 


N=59 


Jnl 


0.0146 


0.0147 


0.0147 


0.0147 


0.0147 


0.0147 


tnl 


0.000308 


0.000312 


0.000312 


0.000311 


0.000311 


0.000311 



Example 2 As a second example, we choose the parameters such that the volume is 
relatively small, V ~ 10^. 

a = 02 = ^, ag = ^, ^2 = 0.04, A3 = 1.2 x 10-^ A4 = 1.2 x 10"^ 



Ai 



-1, A2 = 1,A3 = 0.01, A4 = 0.01,W = 1,^ = 35,7 = 2.65 x 10"^. 



ri(0) = 1781.356, r2(0) = 51.039, r3(0) = 282,T3(0) = 285, 

fi(0) = f2(0) = 0,f3(0) = -1.40948 X 10"^f4(0) = -1.21344 x 10"^. 

The inflation lasts for = 62.5 e-folds. 

One can check that the above explicit results are consistent with the conclusion in 
section (4.2), that Jnl ~ 0{r]). 

^The different values of A'' are computed at correspondingly different, constant, values of the Hubble 
parameter, He- 
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N=20 


N=30 


iV=40 


N=50 


iV=55 


N=59 


Jnl 


0.0171 


0.0167 


0.0148 


0.0108 


0.0121 


0.0170 


tnl 


0.000423 


0.000429 


0.000465 


0.000464 


0.000333 


0.000423 



5. Numerical Methods 

5.1 The Finite Difference Method 

Numerically it is straightforward to solve the equations of motion for the background fields 
without applying slow-roll approximation [10]. The advantage of the numerical method is 
that we do not need to rely on slow-roll approximation (although we still need to assume 
slow-roll at Hubble exit [25]) and no assumption about the shape of the potential is needed. 

We will use the finite difference method [42] to calculate the derivatives of = 
N{(j)i, ...,4)n', Hc)^ up to the second order beyond the slow-roll regime. 

First Order derivative The finite difference method gives 

Ni = ^[N{<Pi, ^^ + h„ ..4n) - N{^i, ct>^ - h„ ...0„)] + 0{h^) (5.1) 
Second Order derivative When i = j, 
Nii = ^[N{(ki, (ki + hi, ..4n) - 2N{cPi, </.„) + N{<l)i, (^i-hi, ..4n)] + 0{h^) (5.2) 
and when i ^ j, 

Nij = ^^^^ [jV((/'i, ■■.,(t>i + hi, ...,(f>j + hj, ...4>n) - N{4>i, ...,4>i + hi, - hj, ...(f>n) 

- N{(pi, (j)i - hi, (pj + hj, ...cpn) + N{(pi, (pi - hi, 4)j - hj, ...</>„)] + 0{h^) 

(5.3) 

Once the NiS and iVjj's are obtained, we are ready to calculate the non-gaussianity using 
the (5 A^- formalism discussed previously^. 

5.2 Example 

We numerically solve the background equations of motion for the model introduced in 
Section 3. Then we use the 5N formalism to calculate the non-linear parameters f^i and 
Tnl- The parameters in the potential (3.15) and initial conditions are chosen the same as 
in Section 4.3, Example 1. 

*Here are understood to be the field values at the Hubble exit. 

^We do not compute the third order derivative of N, and thus Qnl, since the term containing it is 
proportional to 0{h?) which is very small and the error bars can be relatively large. 
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N=20 


N=30 


iV=40 


N=50 


iV=55 


N=59 


iV=60.4 


Inl 


0.00874 


0.0125 


0.0142 


0.0143 


0.0143 


0.0143 


0.0143 


tnl 


0.000300 


0.000274 


0.000292 


0.000295 


0.000296 


0.000362 


0.000346 



These are very close to the results obtained by the analytical method in Section 4.3. 
Remarkably, notice that the values of the nonlinear parameters does not change much near 
{N = 59) and after {N = 60.4) the end of inflation when slow-roll condition breaks down. 
It is reasonable to suspect that the non-gaussianity evolves very slowly through inflation 
and even preheating era. In practice, we may just use the separable potential method to 
compute non-gaussianity under slow-roll condition and use the result as an approximation 
to those in regimes beyond slow-roll. 

6. Curvaton Scenario 

So far we have only considered the inflationary scenario in which we assume that the 
non-gaussianity is generated by the inflaton. However, it is necessary to investigate the 
possibility of a curvaton scenario which does not affect the dynamics during inflation but 
may play a major role in the oscillation stage. 

6.1 Curvaton evolution 

In a multi-field inflationary model, there will in general be several light fields, with one of 
them, called the inflaton, (p, dominating the dynamics of inflation. Other light fields, on 
the other hand, have very little effect during inflation and are usually neglected. However, 
under certain circumstances, a light field other than the infiaton may be identified as 
the curvaton, a, which sometimes generates significant non-gaussianity after the end of 
inflation [43] [44]. 

After the end of inflation, the inflaton quickly starts to oscillate about its potential 
minimum. It then decays into radiation (photon) when its decay rate T(f, ~ H, where 
the decay rate can be calculated once the coupled Lagrangian is given. During the 
oscillation process, if > To-, the inflaton will decay flrst, leaving the curvaton as the 
only light field^. Right after the inflaton decays into radiation, the curvaton energy density 
is still subdominant. However, the massless radiation decreases faster, ~ 4r, than the 
massive particles associated with the curvaton, ~ as the universe expands. Thus the 
relative energy density of the non-relativistic curvaton may increase until it decays into 
radiation, at which point it may even dominate the total energy density. 

6.2 The existence of the curvaton 

We now turn to address the question whether the curvaton scenario can occue in the class 
of models constructed in section 3. For simplicity, we assume that during inflation all the 

®We assume that there are no other light fields, such as those associated with cold dark matter, etc. 
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fields, except (pi and stay close to their VEVs and are thus heavy. We can write the 
potential as 

V^Vo + Vi+Vn (6.1) 

where Vi and Vn are potentials for (pi and (pn, and Vq is the (almost) constant part of the 
potential due to the (almost) frozen (pi, with 1 < i < n. 

If the curvaton exists, its mass must be less than the Hubble parameter. Thus the 
quadratic potential Vi should be small compare to Vq. The inflaton 0„ is displaced far 
away from its VEV, and its potential is suppressed by orders of y 

negligible if compared to Vq. Here /3 depends on how far away (pn has been displaced from 
its minimum. 

As a result, (6.1) is dominated by Vq 

V^Vq = Pq^ (6.3) 



where 

^'o = -.5]aA,(r,)i + -e + 7V (6.4) 



2 

i=2 



(r.) (InV - Inc.) , C, = ^-^^ (6.5) 



where (rj) is the value of the iih. moduli at its minimum and the uplifting parameter 
7V ~ 0(1), see Section 3. 

Near the potential minimum, the masses of the canonicalized fields, (pi and (pi,i > 2, 
given by (3.16) and (3.17) are 

mf = gi^, (6.6) 

^^=Qi^^ ^>2. (6.7) 

where the coefficients Qi are given by 

^1 = -f ^ ^""^^ ' + y ^ + ^^-^^ 

= + 4a-^ (r,) + 4a-2 (r^)' (6.9) 

Note that the mass of the inflaton, after inflation is 0(V~^) while during inflation the 
mass goes like 0(V~^~'^) 
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As expected, the fields (j)i,2 < i < n — 1 are heavier than the Hubble parameter in the 
large volume limit 

^^^-^>H^^l^O-^ (6.10) 

by (3.27). 

Since the field (j)i is our candidate for the curvaton field, it should be lighter than the 
Hubble parameter, i.e., 

< Qi < ^Po, (6.11) 



More explicitly, 



61 " 



si 3 79 17 

J^QAi(ri)2 +-aA„(r„)2 > ^ + (6.12) 



4 

i=2 



^^f;aA,(r,)i <^^ + 67V (6.13) 



4 

i=2 

From (6.12) and (6.13) it then follows that 

3 1 20 

aA„,(T„)2 > -e- -7V (6.14) 

In the simpliest setup where all the fields (pi,i > 2 are identical in parameters(Aj, Oj, 
etc), we get from (6.12) and (6.13) 

n<1.55 (6.15) 

which is not possible since the integer n has to be at least 2. Hence a different setup than 
the simpliest one is needed to satisfy the curvaton condition. 

Equations (6.12) and (6.13) are the necessary conditions for the existence of the cur- 
vaton scenario. They are very restrictive, however, as shown by the simple example above, 
and some fine tuning is required to satisfy (6.12) and (6.13). 

6.3 The decay rate 

It has been shown that the there exists a Lagrangian arising from D7-branes wrapping a 
small four-cycle, whose complexified volume is the modulus r, 

where A^ is the coupling of r to the gauge field (photon) [45] . 

The other parts of the Lagrangian can obtained by quadratic expansion around the 
potential minimum. By canonically normalizing the kinetic terms and diagonalizing the 
mass matrix terms 

Co = -V^in + Idf.i^.d'^ij' - ^mjilj.i;' (6.17) 

where 'i/'i are the canonically normalized fields and also eigenfunctions of the mass matrix. 

In what follows, we consider a model consisting of multiple moduli: ti,T2, t„, where 
Ti is the large four-cycle and all other moduli are small. Typically, most of the small cycles 
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are close to their vevs and thus are heavy during inflation. Only the inflaton is displaced far 
from its vev. Let the inflaton be Tn- So the relevent moduli here will be ri and r„ and other 
moduli play the role of stabilizing the potential. Starting from the Lagrangian in (3.13), 
it is possible to simultaneously diagonalize the kinetic terms and the mass matrix terms 
under the assumption that the mass matrix is independent of the fields, which is a good 
approximation close to the minimum of the potential. For simplicity, let us first diagonalize 
the kinetic terms using the field transformations (3.16) and (3.17). The Lagrangian reads 

1- - 1" - 1 -2l -2 
jCo = - Vmin + -d^4>ld'^(f>l + -d^(f)nd'^(j)n - ^mjcpi - -m?n(j)n - mf„0i(/)„ 

where 

- o , r., - Ql,n—5- (6.1«j 
0(l)ld(j)n V2 

Qi,n ^ -?>\plaXaan (r^)^ , (6.19) 

and the effective fields (j)i = (j)i — {(pi) represent the oscillation amplitude of the field about 
its potential minimum. 

From now on we omit the hat over (p. Following [46] , we calculate the eigenvalues and 
eigenvectors of the mass metrix 

^2^/m|m2\ ^^^^^^ 



M'^Vi = eiVi, i = l,2. (6.21) 



where the eigenvectors Vi are normalized such that vfvj = 6ij. 
The transformation takes the form 



j = (^'i) V'l + (^^2) 4^2 (6.22) 



In the large volume limit, we find that 

~ + (V~^) ^ V"! (6.23) 
(t)2^0 V'l + ^^2 ^ 4^2 (6.24) 

which means that we can approximate the Lagrangian by 

Co ~ -y„,in + ^df,<Pid''(l)i + ^df^^nd^^cPn " ^mf^ - ]^ml4>l (6.25) 
Similarly, we can rewrite the Lagrangian for the gauge sector [45] 



-18- 



where k is & normalization factor. We set t = Tn since the D7 branes only wrap the small 
four-cycle [41]. In terms of the (approximately) canonically normalized fields (pi, the 
gauge field (radiation) Lagrangian (6.26) takes the form 

1 " A 



1=1 



which corresponds to 



k={Tn)-\ (6.28) 



2 ' ^" V4aA„ 
The complete Lagrangian reads 



1 



\/6 / 3V \2 3 

A0i = ^, A^„= -— (r„)-4, A^, «0, i = 2,...,n-l. (6.29) 



From the Lagrangian (6.30), it is straightforward to get the decay rates 

i.e., 

^1 1 



3 



(6.31) 



r0i-5.77 — ioo_°.r2 ^1,.9 (6.32) 



Thus, 

Ti ~ At « r„ ~ ^ (6.34) 

This indicates that the curvaton (pi decay indeed occurs some time after decay of the 
inflaton. The amount of non-gaussianity generated by the curvaton is determined by its 
relative energy at the time of decay, which will be shown in the next section. 

6.4 The nonlinear parameter 

The curvaton starts to decay when the Hubble parameter drops below the decay rate of 
the curvaton ^ 

H ~ Tfj^^^ (6.35) 



^In what follows, the curvaton in our scenario, cj>iy is relabeled a to conform with previous work on 
curvatons in the literature. 
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Using the sudden decay approximation (assuming the decay happens instantaneously), 
the nonhnear parameter of the curvaton perturbation can be shown to be [24] [47] 



^ _ 5 5rde 

4rdec V ff'V 3 



fNL = j^{l + '^]-^--^ (6.36) 



where the dimensionless ratio 



^dec = ^- —J- (6.37) 

and Pcrdcc and p^dcc are the energy density for the curvaton and radiation at the time when 
the curvaton decays, respectively. 

The function g characterizes the dependence of the curvaton, cj(= (pi), at the beginning 
of its oscilation, on its value at Hubble crossing, a*, i.e., a = g{a^) . Assuming the absence 
of the nonlinear evolution of the curvaton, we have g" = and 

fNL = ^-\-^-^ (6.38) 
4rdec 3 6 

If the curvaton energy density dominates the total energy density when it decays, the 
corresponding nonlinear parameter is 

fNL--\ (6.39) 

On the other hand, if rdec ^ then f^L ^ 1- Note that if rdcc ~ 0.58, Jnl ~ 0. 
The initial energy density of the curvaton cr, when it begins to oscillate, is 

Po - \mla^ (6.40) 

where a is the oscillation amplitude of the curvaton. To estimate its value, we use the 
arguments similar to [41] [48]. Analogous to the Hawking radiation in black holes, the 
quantum fluctuation 5a of the light field a during inflation in de Sitter space has the power 
spectrum [49] 

|2\ / \ rr2 



= (IHO = ^i^J =Tfi (6.41) 
where Th is the Hawking temperature and the label -k denotes the Hubble exit; and 

Hi ^ \P.^ (6.42) 
by (6.3). It indicates that the amplitude of quantum fluctuation 

5a^TH = ^ (6.43) 

The amount of quantum fluctuation is comparable to the classical (slow roll) motion 
when 

5a ~ a5U ~ i^H'^ (6.44) 
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where the slow-roll condition of the light field a has been used, and 61^, = is the 

change in time during one e-fold. We view the onset of the quantum regime as the time 
when the oscillation takes place. The typical (initial) value of a constraint by the quantum 
fluctuations thus satisfies the conditon 

dV n 

^ ^ Ht (6.45) 
Since the potential is quadratic under assumption, the value of a reads 

3 

a^Ar (6.46) 

where K = Vb ~ Sif^. 

The initial ratio between the curvaton energy density and the total energy density is 



where 

and we set Mpi = 1 as usual. 
Since 

we have 



i^in = (6.48) 



ml = ml = Qi^, (6.49) 

Now assume that the oscillation stage only lasts for a few e- folds (AA^ = J^'^ Hdt). 
We should have <C e~^^ < 1. Then, under the sudden decay approximation, the ratio 
rdec of Eq. (6.37) can be related to J7in by [50] 

Tdec ^ |^^m(l - f^in)^^ {~Y^ 

7 i 'ii 



(6.51) 



4 V r 



67rPo3vyg2 1 

5 3~ 
2 V2 



In terms of iJ*, 



Tdec^Wi^ (6.52) 

where we temporarily restore the Planck mass. Notice that the Hubble parameter during 
inflation is generally much smaller than the Planck mass. Given the fact that Wq should 
not be very large, the ratio r^ec can be quite small. By (6.38), this will give rise to a large 
positive /nl > 1, see also [41]. 
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7. Conclusion 



In the first half of this paper, we discussed the inflaton scenario. Both the analytical 
method (separable potential) and the numerical method (finite difference) have been used 
to calculate the nonlinear parameters via the 6N formalism. These two methods agree with 
each other very well in the string inflation model introduced in Section 3. We conclude 
that, although the analytical method is only valid under slow-roll, it can be a good estimate 
for regimes beyond slow-roll. The nonlinear parameters we get are very small which is 
typical for slow-roll inflationary models. The results do not vary much for different volume 
sizes. Instead, as shown in Section 4.2, they mostly depend on the slow-roll parameters. 
The amount of non-gaussinaity {/nl ~ 0.01) are unlikely to be detected by the major 
experiments that will be carried out [51] [52]. However, a certain type of cosmological 
probe has been proposed [53], by which this amount of non-gassianity is reachable. 

After discussing the inflaton scenario, we present the conditions for existence of the 
curvaton scenario. To satisfy those very restrictive conditions, the parameters need to be 
fine tuned appropriately. The non-gaussianity in the curvaton scenario is given in terms 
of the parameters in the potential and the Calabi-Yau volume. A rough estimate shows 
that the non-gaussianity can be quite large for a typical setup. These large non-gaussianity 
effects in curvaton scenarios have also recently been computed in closely related type IIB 
flux compactifications [41]. 
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